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ABSTRACT 

We study the strong gravitational lens statistics of triaxial cold dark matter (CDM) halos occu- 
pied by central early-type galaxies. We calculate the image separation distribution for double, 
cusp and quad configurations. The ratios of image multiplicities at large separations are con- 
sistent with the triaxial NFW model, and at small separations are consistent with the singular 
isothermal ellipsoid (SIE) model. At all separations, the total lensing probability is enhanced 
by adiabatic contraction. If no adiabatic contraction is assumed, naked cusp configurations 
become dominant at » 2.5", which is inconsistent with the data. We also show that at small- 
to-moderate separations (< 5") the image multiplicities depend sensitively on the alignment 
of the shapes of the luminous and dark matter projected density profiles. In constrast to other 
properties that affect these ratios, the degree of alignment does not have a significant effect on 
the total lensing probability. These correlations may therefore be constrained by comparing 
the theoretical image separation distribution to a sufficiently large lens sample from future 
wide and deep sky surveys such as Pan-Starrs, LSST and JDEM. Understanding the correla- 
tions in the shapes of galaxies and their dark matter halo is important for future weak lensing 
surveys. 



1 INTRODUCTION 



Statistics of strongly lensed images promise a wealth of informa- 
tion about luminous and dark matter distributions and the corre- 
lations between them. For example, these statistics have yielded 
gen eral constraint s on th e radial mass profil e of ea rl y-type galax- 
ies jMaoz&Rixl Jl993h. iKeeton & Madaul JioOlh. iwvithe et al l 
1 2001 ), iKeetonl feOOlbl). [Rusin & Mai feOOlh. iTakahashi & Chibal 
1200111, ISarbu et alJ feOOll). lOguril j2002h. lOguri et alJ feOoUT 
iHuterer & Mai fe004h T Tkuhlen et al.l fe004h ). Lensing by non- 
spherical mass distributions shows rich phenomenology that may 
be used to probe the distribution of dark matter in galaxies and 
clusters. The shape and concentration of the lenses Mass distribu- 
tions that are approximately spherical produce a higher fraction of 
double-image systems, whereas more ellipitical distributions pro- 
duce more quadruple-image systems. In rare cases, extreme ellip- 
ticities can produce strongly magnified three-image configurations 
known as "naked cusp" systems (or simply "cusps"). Studying the 
statistics of different image multiplicities (quad, double, and cusp 
lenses) could help us constrain the shapes of dark matter halos, the 
environment of and substructure in these halos and as we will show 
the correlations between the shapes of the luminous and dark matter 
in galaxies. 

The correlations in the shapes of the galaxy and the dark mat- 
ter halo it resides in is important for four reasons. (1) As we show, 
these correlations are important to include to make accurate pre- 
dictions for the strong lensing probabilities. (2) A large sample of 
strongly lensed images will provide a rich data set that can help 
us understand galaxy formation better. One of the most interesting 
aspects of galaxy formation, and perhaps the least well studied, is 
the correlation in the shapes of the luminous and dark matter com- 



ponents. Lensed image statistics can provide constraints on this 
aspect. (3) An important systematic effect for weak lensing sur- 
veys is the spatial correlation in the shapes of the galaxy shapes, 
which arise is if galaxy sha pes are correlated with t h eir dark matter 
halo a n d halo environment ([Croft & Metzlej 1 2000h. Heavens et al.l 
(2000), Cat elan et all feOO lLlCrittenden et alJfeOOlh . lJinsI d2002l) . 
Mandelbau m et alJ ( feOOfjh . iHirata et all j2007h ). All future weak 
lensing surveys will cover a large fraction of the sky at great depth 
and also reveal a great number of strongly lensed quasar and galaxy 
images. With an accurate framework in place, one may use this 
strong lensing data set to understand systematic effects stemming 
from luminous and dark matter correlations for weak lensing bet- 
ter. (4) Strong lensing systems are a lso important for un derstanding 
fundamental physics. As shown bv lBolton et alJ | |2006|) . these sys- 
tems may be used to constrain the post-newtonian gravity param- 
eter 7. Constraining this parameter also directly contrains theories 
of modified gravity since in this context 7 i s the ratio of t he tw o 
scalar potentials i n the perturbed FRW metric teertschingeJ ( l200fj) : 
IZhang et al.U2007h ). 

The distribution of image multiplicities for elliptical lenses 
has been modeled in two regimes: small image separations (< 1"), 
in which galaxies become the dominant lensing mass, and large 
(cluster-siz ed) separations (> 10"), in which dark matter becomes 
dominant. iRusin & Tegmarld ( 1200 ll) modeled galaxy lenses with 
a singular isothermal ellipsoid (SIE) mass distribution, compar- 
ing to lenses from the Cosmic Lens All-Sky Survey (CLASS) and 
lodrell-VLA Astrometric Survey (IVAS). They showed that the ob- 
served number ratio of quads-to-doubles is greatly underestimated 
by the predictions of the SIE model (dubbed the "quad problem"). 
While this problem has not been solved rigorously, it has been sug- 
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gested that the effect of massive substructures could sufficiently 
increase the expected quad fracti on to be consistent with obser- 
vations dCohn & Kochanekl d2004h l. Other effects which boost the 
lens probability for quads inclu de misalignment of the centers of 
galaxies with their host halos ( Ouadri et al j j2003[ ) ) and external 
shear from lens galaxy environments (Hute rer et al.1 d2005l) ). 

Large separation l enses, on the other hand, were modeled by 
lOguri & Keetorj using the triaxial NFW model based on 

simulations by Ijing & Sutol (2002). With an inner slope a = l.( 
(pure NFW), the image multiplicities are shown to be remarkably 
different from those of galaxy lensing: cusp systems have the high- 
est probability, followed by quads, with doubles being least com- 
mon. However, because the overall lensing probabilities are orders 
of magnitude smaller than that of galaxy lensing, such dark-matter 
dominated lenses are relatively rare — the only observed lens in this 
regime is a quad a t 14.62", found in the Sloan Digital Sky Survey 
dlnada et al.1 d2003h ). 

In contrast to the paucity of confirmed large-separation lenses, 
a significant fraction of observed gravitational lenses lie in the in- 
termediate region between these two extremes, with image separa- 
tions roughly between 1" and 10". The total lensing probability in 
this regi on has been s tudied for t he case of spherical mass distri- 
butions dOguril d2006h , lM3 d2003f ), lKochanek & White] d200ll) ). but 
the distribution of image multiplicities for nonspherical lenses has 
not yet been derived. We note that none of the single-component 
elliptical models mentioned thus far can be applied in this regime, 
since both the dark matter and the central galaxy contribute signifi- 
cantly at these separations. Furthermore, any lensing model includ- 
ing baryons and dark matter must take into account the modifica- 
tion of dark matter halos due to baryon cooli ng, which increases 
the co n centration of th e inner parts of halos dKochanek & Whitel 
d200ll) . lKeetor]d2001bl)). Althou gh a multi-component model was 
employed by Mo ller et alj d2007l) using galaxies with ellipticity, the 
dark matter was modeled by spherical NFW halos, ignoring adia- 
batic contraction and the triaxiality of dark matter halos. 

In this paper we derive a distribution of image multiplicities 
that will be approximately valid for all observed image separations, 
including the intermediate region. However, to accurately model 
this region, we must combine the triaxial NFW model, adiabatic 
contraction, and a central lens galaxy. The lensing properties then 
become sensitive to the ellipticities, concentrations, and orienta- 
tions of both galaxy and halo, in addition to the fraction of the total 
mass contained in the galaxy. This requires a more complicated 
calculation of lensing probabilities than those in the literature. De- 
spite the complexity of the models considered in this work, several 
important factors (considered prevoiously in isolation) have to be 
neglected. In particular, we highlight the absence of dark matter 
substructure in our model and the effect of the lens galaxy environ- 
ment. We will later discuss the expected impact of these factors. 

We find several interesting features that suggest new methods 
of constraining lens mass distributions. For example, the ratio of 
image multiplicities is found to be sensitive to the alignment of 
the major axes of galaxies with that of their host halos. Further, 
the transition from small- to large-separation lenses is marked by 
a "crossing" of lens probabilities when the expected lensing rate 
of quads overtakes that of doubles, and likewise for cusps. Under- 
standing the precise location of these turnover points may provide 
further means of constraining galaxy-halo mass distributions once 
the sample of observed lens systems becomes sufficiently large. 

The paper is organized as follows. In Sj2]we will present a 
general form for the image separation distribution, and discuss the 
relevant model parameters. Section [3] will describe the lens model 



we used and outline our approach toward applying the adiabatic 
contraction model to triaxial halos. Section [4] covers the details of 
the calculation, and results are given in §5\ Finally we summarize 
and discuss our results in ^6] 



2 DISTRIBUTION OF IMAGE SEPARATIONS 

The differential image separation distribution is usually defined as 
the probability that a source object at redshift z s wil l be multiply 
image d with an image separati on 0. This is written as dTurner et al.1 



( 1984). ISchneideretal.l 



iparatioi 
(2006)) 



dP C ,cdt , C dn dM 

— = J dz(l+zf—D- a „ g (z) J dXiP(Xi) B<r( & , Xi ,z)-- 

M(@, Xi ,z) 

Here Bcr is the biased lensing cross section (explained in 5j4]l in 
units of solid angle, and n is the comoving number density of halos. 
The lensing optical depth is given by the integral of this distribution 
over all image separations; in this paper we will simply refer to it 
as the total lens probability, though always with respect to a given 
source redshift z s - The Xi Me the parameters relevant to lensing, 
whose distributions p(xi) are integrated over, and which depend on 
the particular lens model chosen. In this study we take our model to 
be an adiabatically-contracted NFW dark matter halo with a central 
early-type galaxy represented by a Hernquist profile. The lensing 
parameters are then 



X = (c e ,a/c, b/c, 9, (f>, f b , x b , q H , 9 q ) , 



(2) 



where c e is the triaxial concentration parameter (defined below) of 
the original dark matter halo before adiabatic contraction, ajc, b/c 
are the axis ratios of the halo, (9,<f>) are the orientation angles of 
the halo, f b is the baryon mass fraction, x b = r b /r s gives the scale 
radius of the Hernquist profile, q H is the projected axis ratio of the 
Hernquist profile, and 9 q is the angle between the projected major 
axis of the Hernquist profile and the projected major axis of the 
dark matter halo. 

The total mass of the lens was not included among the pa- 
rameters listed above because for a given image separation 0, lens 
redshift z and set of lensing parameters Xu me mass is uniquely 
fixed. 



3 THE LENS MODEL 

In this section we describe the model for the lenses (elliptical galax- 
ies and clusters). We split the discussion into four parts dealing 
with the dark matter halo, the baryonic component, effect of baryon 
cooling on the dark matter halo and finally correlations between the 
luminous and dark matter mass distributions. 



3.1 Triaxial dark matter halo 

Throughout this paper we assume a concordance cosmology with 
Q. M = .3, n A = 0.7, h = 0.71 , and cr 8 = 0.9 and use the Seth- 
Tormen dSheth&Tormer]d2002l) ) mass function of dark matter ha- 
los. 

We assume the de nsity profile of a triaxial dark matter halo to 
be dJing& Sutol d2002h ). 



p(R) 



(R/R )(l+R/R ) 2 



(3) 



where 
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<x 2 



b 2 



(a < b < c). 



(4) 



The triaxial concentration parameter is defined in ljing & Sutol 

(2002) as c e = R e /Ro, where R e is defined so that the mean density 
within an ellipsoid of major axis radius R e is A e Q M (z)p C nt(z) with 
A f = 5A rjr (c 2 /ab) 0J5 . The characteristic density is then written in 
terms of the concentration parameter as 



See = 



A e £l M (z)c 3 e 



(5) 



3 [ln(l + c e ) - cj(l + c e )] 

We use fitting functions from ljing & Sutol {2002) for the dis- 
tribution of axis ratios and concentrations (nearly iden tical distribu- 
tions f or the axis ratios were found in simulations bv lAllgood et al.l 
(2006)). The angular orientations are of course entirely random 
(p(&) = (sin6>)/2, p(<t>) = l/ 2?r). The fitting func tion for the median 
concentration parameter in Ijing & Sutol 1 20021) has a dependence 
on the axis ratio a/c; however the scatter around their fit is quite 
large, and in practice this only affects the image separation distri- 
bution significantly at large separations (© > 30"). Therefore we 
will neglect the a/c dependence in this paper. 



3.2 Baryonic component 

For modeling the projected surface density of the galaxy we use an 
elliptical Hernquist profile defined by 



1H 



A' 2 



1H 



(6) 



where £ //(r) is the projected circular Hernquist profile (see lKeetonl 
(2001a) for the analytical formula) and q H is the projected axis ra- 
tio. The ellipse parameter £ is normalized so that the mass enclosed 
inside a contour of constant surface density is independent of q H \ 
this ensures that, when averaged over all angles, the circular pro- 
file S H (r) is recovered. This method of normalization is only ap- 
proximately correct, since the normalizatio n changes if there is a 
tendency toward oblate or prolate galaxies (Chae (2003)). To ob- 
tain an average projected profile of axis ratio q, one must average 
over the expected distribution of axis ratios and viewing angles that 
produce the projected axis ratio q. However, our approximation is 

[ ood if we a s sume equal number s of oblate and prolate galaxies 
Chae] j2003h : lHuterer et alj J2005h ). 

In order to relate the baryon fraction and scale radius to 
the mass and concentration of the host halos, we adopt the 
mass-luminosity relation proposed by IVale & OstrikeJ J2004) in 
terms of bj-band luminosity. We convert luminosities to stellar 
masses by assuming a universal stellar mass-to-light ratio of T = 
3.0hjoM Q /L Q . We use the scaling relation between g-band luminos- 
ity and effe ctive radius R derive d from 9,000 early-type galaxies 
in SDSS bv lBernardi etall J2003I) : 



(7) 



10" 



10°- 



2 h 7 2 kpc 



For a spherical halo, the effective radius is related to the scale 
radius of the Hernquist profile by r b = 0.55 IRq. However, this rela- 
tion changes for an elliptical profile since R refers to the observed 
circular effective radius. Therefore we use a relation of the form 
r b = 0.551 f(q H )R , wh ere f(q ) is de termined by numerical in- 
tegration of X w (£). As in lOguril feOOd) . we neglect the difference 
between bj-band luminosity and g-band luminosity for simplicity. 



3.3 Baryon cooling 

Adiabatic cooling of baryon leads to a increase in the concentra- 
tion of dark matter halos. We incorporate the effects of baryon 
cooling on the dark matter halo by adopting a modified adiabatic 
co ntraction (AC) model. Th e model is based on the original idea 
of Hkimenthal et~ai] dl986l) . which assumes spherical symmetry 
and uses conservation of angular momentum to calculate the re- 
sponse of dark matt er to baryon infall. The modified AC model of 
iGnedin et al] J2004I) incorporates the eccentric orbits of particles, 
though spherical symmetry of the overall density profile is still as- 
sumed. From simulations based on this model, they presented a 
series of fitting functions that map the radius of a particle after con- 
traction to its radius before contraction. The final density profile 
as a function of radius can then be found by inverting the map- 
pi ng via interpolation . The only caveat here is that the simulations 
in lGnedinetal.l ( l2004l) apply only to spherical NFW halos, whereas 
we are assuming a triaxial model. To describe the density profile in- 
dependent of the effect of AC on halo shape, we apply the modified 
AC model with the replacement r — > R, where the triaxial radius R 
is defined as above. We expect this to be a good approximation as 
long as the axis ratios are not too small. 

In addition to making halos more concentrated, adiabatic 
contraction is also expected to make them more spherical, driv- 
ing the a xis ratios closer to one especially in the inner parts of 
the halo dDubinskil jl994h ). This effec t is seen in s imulations by 
iKazantzidis et al.ld2004l) (see also lKazantzidis et alj J2006I) ). How- 
ever, the impact of these results on lensing is not well understood. 
The image separation and cross sections are particularly sensitive 
to the density profil e in the inner parts of t he halo. However, the 
CDM simulation of Ijing & Sutol J2002h and I Allgood et all 120061) 
actually show that, in the absence of baryons, the axis ratios tend to 
decrease slightly toward the center of the halo. Whether the inner 
axis ratios of the dark matter will end up greater than, less than, or 
roughly equal to the initial average overall axis ratios is unknown. 
Therefore, for most of the paper we will assume the distribution 
of axis ratios after contrac t ion ar e not substantially different from 
those given in |jing& Sutol J2002h . Later we will investigate the ef- 
fect on lensing if the axis ratios are increased following adiabatic 
contraction. 



3.4 Correlations between the luminous and dark matter 

Luminous matter cools and forms stars in the gravitational poten- 
tial well of the dark matter and hence it is reasonable to assume 
that the shapes of the dark and luminous matter mass distributions 
should be correlated. There theory explaining this correlation has 
not been worked out. Here we ask how these correlations might 
affect the lensing probabilities. In particular, we correlate the ellip- 
ticity and projected major axis of galaxy to that of the dark matter 
halo it resides in. Previous studies h ave suggested a possible hig h 
degree of alignment for disk galaxies jvan den Bosch et alj j2002)). 
However, there is not a clear consensus. We will study the image 
separation distribution for the extreme cases, in which the ellip- 
ticities are completely correlated vs. uncorrelated, and completely 
aligned vs. misaligned. We note that the case where the there is no 
correlation in the shapes of the dark and luminous components will 
lie in between the two extreme cases we consider. For the case of 
completely uncorrelated ellipticities, for simplicity, we will fix q H , 
the ellipticity of the galaxy (Hernquist profile), to a particular value 
and integrate over the distribution of dark matter halo ellipticities. 
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When making predictions for a given survey, our results have to be 
generalized by averaging over q H using the observed distribution. 



4 CALCULATING THE IMAGE SEPARATION 
DISTRIBUTION 

To do the lensing calculations, for each lens we define dimen- 
sionless coordinates in the lens plane i = ^(z,z s )/^o(z,z s ), where 
£ has units of distance in the lens plane and £o(z, z s ) is taken to 
be the Einstein radius of the corresponding spherically symmet- 
ric lens at a redshift z and source redshift z s - (The angular separa- 
tion © is given by f = &D L (z), where D L is an angular-diameter 
distance.) We also define corresponding coordinates in the source 
plane, X = ff(z,z s )/r]o(z,z s ) where rj = g D s /D L . The cross sec- 
tions and image separations in these units will be denoted by & and 
0, respectively. 



4.1 Projected mass density profile 

Before solving the lens equation we must first obtain the projected 
mass density profile, Yif). This is generated numerically by taking 
the triaxial adiabatically-contracted NFW profile mentioned above, 
assuming some orientation angles 6 and <p, and integrating along the 
line of sight. However, this method can be simplified considerably 
by a judicious change-of-variables. We give a brief su mmary of the 
metho d her e; for more d etails we refer the reader to lOguri et al.l 
(2003) and JOgurill2004 pg. 40), where the method is applied to 
triaxial NFW halos. The triaxial radius R (defined above) is written 
in terms of the observer's coordinates (x',y',z r ), with z' axis along 
the line of sight, as 



V/z' 2 + g(x',y')z' +h{x',y') = -yjz', 2 + f 2 



(8) 



where /, g, h are all functions of the axis ratios and orientation an- 
gles, but / has no dependence on x',/. Further z', = -\/7(z' +g/2f), 
f 2 = h - g 2 /4f and with this change of variables, we have 



(9) 



The ( parameter can be shown to correspond to the pro- 
jected isodensity curves, the two-dimensional analogue of R de- 
fined above. It can be written as 



,2 



, x y 
? = — + — 

q 2 % 



(10) 



where the x-axis is taken to lie along the projected major axis of 
the isodensity curve, and q x and q y are complicated functions of 
the axis ratios and orientation angles. We can see from the equa- 
tions above that the result of changing the axis ratios and orien- 
tations is to scale the projected density by a factor F z = 1/ ^ff, 
and scale the x- and y-coordinates by the factors q x and q y respec- 
tively. Therefore, we need only calculate the above integral for the 
spherically-symmetric case, and then apply these scalings to obtain 
the projected density profile for any combination of axis ratios and 
orientation angles. This vastly simplifies the calculation. Finally, 
we superimpose this with the Hernquist profile described above, 
and divide by the critical lensing density Y. cril to arrive at a dimen- 
sionless "kappa" profile k(x). 



4.2 Lensing cross-sections 

After generating the kappa profile, we then calculate the specific 
cross sections for lensed sources of different image multiplicities, 
depending on whether they are double, quad, or naked cusp image 
configurations. The cross sections and image separations are cal- 
culated by a Monte Carlo method where a thousand sources are 
placed randomly within the caustics of the gravitational lens. For 
each source point we solv e the lens equat ion numerically, using 
the integration formulas o f ISchramml (fl99Ch for the deflection and 
magnification of elliptical mass profiles to o btain the corresponding 
images and magnifications (Keetonl( [200T3) ). To find the images we 
use an adaptive radial grid where the grid cells recursively divide 
themselves several times, with two more subdivisions when a cell 
is in the vicinity of a critical curve. Images are then found by map- 
ping each trapezoidal cell to the source plane and testing whether 
the source point falls inside the mapped region. (Each cell is actu- 
ally split into triangles first to avoi d the ambiguity of having convex 
vs. concave cells; see lBartelmannl d2003l) for details.) 

We exclude all lens systems with a flux ratio between the 
brightest two images greater than 15:1, in accord with the flux ra- 
tio selection function of radio searches using the VLA. In practice, 
such large flux ratios are common only in the isothermal case with 
double images. For the adiabatically contracted model, less than 
0.1% of lensing events result in such large flux ratios, whereas the 
fraction for an isothermal lens is roughly 15-20%. 

The image separations, defined here as the maximum separa- 
tion between any two pairs of images in an image configuration, are 
then averaged over all sources to obtain the average image separa- 
tion 0. The cross sections can be obtained by taking the fraction of 
sources which produce a given image configuration (quad, double, 
or cusp), multiplied by the total cross-section. In addition, if we 
also weigh the sources by their magnifications appropriately, we 
can compute the magnification bias. Assuming the sources follow 
a power-law luminosity function 0l(L) oc L~P (generally true for 
luminosities relevant to lens searches), the biased cross-section is 
written as 



B& 



I 



d 2 X 



f 



d 2 X^ 



(11) 



where fi is the magnification, and we are integrating over the 
multiply-imaged region of the source plane in dimensionless units. 

There is, however, an ambiguity in this formula: should the 
lens systems be weighted by the total magnification of all the im- 
ages, or by one of the brighter images? The answer depends on 
the particular survey and the method of selecting lens candidates. 
For most surveys to date, such as JVAS/CLASS, the number of ob- 
served quasars is sufficiently small that each quasar can be returned 
to and observed in great detail. In such cases, for small-to-medium 
image separations the brightness of the lens system is seen as the 
total brightness of a ll the images, so the total magnificatio n shou ld 
be used in the bias ( Takahashi & Chiba l200ll). ICen et al. I fl994h ). 
However, for large angular separations, a lens system will only be 
identified if a second image is independently resolved in the sur- 
vey, and therefore the magnification of the second-brighest image 
should be used. 

Looking ahead to future large-scale synoptic surveys using 
telescopes such as LSST, the number of observed lenses may be 
orders of magnitude larger and each lensed object may not be ob- 
served in such detail. Lensed quasars, for example, can be recog- 
nized using image subtraction m ethods that obviate the need for 
extensive follow-up observations dKochanek et alj d2006h ). In that 
case the second-brightest image bias may be more appropriate as 
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the more conservative bias. This will almost certainly be the case 
when searching for galaxy lenses in large surveys such as the Su- 
perNova Acceleration Probe (SNAP). We will consider both types 
of bias in this study. 

From the method above we can obtain the biased cross sec- 
tions for each image configuration type separately as well as the 
total biased cross section. These will be written as <r, where i = 
2, 3, 4, for doubles, cusps, quads, and the total cross section re- 
spectively. Throughout this paper we will assume a power-law 
index /3 = 2.1 for the luminosity function, which is consistent 
with flat-spectrum r adio sources seen by surv eys such as CLASS 
dMvers et alj j2003h . lRusin & Tegmarkl 1200 lh ). 

The principle difficulty in doing the above calculation is that 
the cross sections have a complicated dependence on almost all 
the lensing parameters Xu m addition to the redshift z. However, 
a few simplifications can be made. We can reduce the number of 
parameters by one if we transform the orientation and ellipticity 
parameters to their projected counterparts: 



(a/c,b/c,e,<f>)^(q x ,q,F z ) 



(12) 



Here q x and F z are the same parameters introduced earlier, and 
q = q y lq x is the projected axis ratio. The parameter q x is the ra- 
tio of the projected major axis to the triaxial major axis of the 
halo. It might be objected that we cannot perform this reduction 
since the characteristic density S ce has an explicit dependence on 
the axis ratios a/c, b/c through the triaxial overdensity A e (see 
eq. |3). However, this dependence simply scales the kappa profile 
and thus can be accounted for by redefining the F- parameter as 
F z = {c 1 lab) ' 15 1 -Jf, where / is defined in eq.(8] 

For simplicity we neglect scatter in the mass-luminosity rela- 
tion and halo concentration. We will investigate later (see ^5} how 
the scatter might affect the calculation. 



4.3 Image separations 

Another simplification can be made by first noting that the average 
image separation, in the absence of ellipticity, is very well approx- 
imated by x 2R E where R E is the Einstein radius; or, in the 
units defined above, © ss 2. Introducing ellipticity changes the pro- 
portionality by a certain amount; however, c e , q x , F, and z simply 
scale the Einstein radius without significantly altering this relation. 
We find that the image separations resulting from the dark matter 
halo is fit well by the form, 



® dm (q) = a(M) -J3(M) [expKl - q)y(M)] - 1] 



(13) 



where M is the mass of the halo. Likewise, the Hernquist compo- 
nent is well fit by 



® H (q H ) = a(M) 0.5(1 + , 



(14) 



Both models are conceived so that a(M) = 2 when q = 1. Finally, 
we model the total image separation by a fitting function of the 
form 



= a(M) 



(RE, dm ® dm (q dnl )f m + (R E ,H&n(qH)f (t 



(RE, dm ® d ,„(q d ,„ = m M) + (Re,h®h(«h = nf m 



• i .1/ 

/JIM) 



■(15) 



The advantage of this approach is that, for a given image sep- 
aration 0, we can easily find the total mass of the lens as a func- 
tion of the other parameters, M(c,q,q H ,q x ,F z ), and likewise for 
dM/d®. This reduces the number of parameters by one. Thus we 
are left with the following independent variables on which the cross 
sections depend: 



<r, frjiq, q H ,q x , F z ,z). (16) 



4.4 Projected ellipticities 

To simplify the analysis, rather than taking qu as an independent 
variable we will investigate two cases: first, we assume the ellip- 
ticities of galaxy and halo are strictly correlated; second, we as- 
sume different fixed values of q H . Our motivation is that since 
the ellipticities of the adiabatically-contracted halos are still not 
well understood, the extent to which the ellipticities are corre- 
lated is unknown. Undoubtedly the true distribution of elliptici- 
ties will lie somewhere between these two extreme cases. If we 
assume that q and q H are completely uncorrelated, we can simplify 
the analysis by fixing q H at various values including its mean ob- 
served value. Several studies have analyzed the isophotal shapes 
of E/SO galaxies in various surveys, the distribution of which is 
typically well fit by a Gaussi an: dn/de oc exp[-( e - eo) 2 /2Ae 2 ], 
where e = (1 - q 2 H )l(l + a 2 :l ). iKeeton et alj <fl997h studied galax- 
ies in the Coma cluster and found eo = 0.26, Ae = 0.33, whereas 
lHao et all j2006h studied 847 E/S0 galaxies in SDSS and found 
6 = 0.22, Ae = 0. 14. The Coma distribution has a mean value at 
(q H ) = 0.67, which we will use for the uncorrelated case unless 
otherwise noted. On the other hand, if we make the simplifying as- 
sumption that q and q H are strictly correlated, we can relate them 
directly by q H = q + Aq where Aq is a fixed increment (or simply 
qu = 1.0 if q > 1.0 - Aq). We will investigate both the correlated 
and uncorrelated cases, taking Aq = OA for the correlated case. 

We will calculate dP/d® at specific values of z; these vari- 
ables will be integrated over at the end. To do the integral we first 
calculate the values of a over a regular grid in the remaining three 
parameters (q,q x ,F z ) and then interpolate to find its value at any 
given point. The interpolation is done by an N-dimensional ana- 
logue of bicubic interpolation, whereby the function values, its par- 
tial derivatives, and higher-order mixed partial derivatives are tabu- 
lated on a regular grid (the derivatives being calculated in this case 
by finite differencing across adjacent grid points). The interpolating 
function is a cubic polynomial in the N variables whose coefficients 
are uniquely determined by matching its values and derivatives to 
the tabulated values on the grid. Thus the resulting function is con- 
strained to have the tabulated values on the grid, and also to vary 
smoothly from point t o point. (For a detailed deri vation of the three- 
dimensional case, see iLekien & MarsderJ d2005h .) With cr in hand, 
we integrate over the axis ratios and line-of-sight angles (which are 
the remaining lens parameters in eq.[2} using the Vegas algorithm, 
and finally integrate over redshift to obtain the image separation 
distribution. 



5 RESULTS 

Throughout this paper we plot the image separation distribution for 
a source placed at a redshift of z s = 2.0; the results are qualita- 
tively the same for other source redshifts. First we assume the ma- 
jor axes of galaxy and halo are aligned, with correlated ellipticities. 
Using the "second brightest image" magnification bias, the results 
are shown in fig. [T] The lensing probability for quads overtakes 
that of doubles at 2", whereas cusps become dominant at 5". Note 
that here we have not taken into account the angular resolution of 
surveys, which cuts off the distribution sharply at small image sep- 
arations. The lensing probability for CLASS, for example, cuts off 
at 0.3". 
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Figure 1. Distribution of image separations for double, quad, and cusp im- 
age configurations. Galaxy and dark matter halo shapes are aligned, with 
correlated ellipticities (qu = q + 0.1). The "second brightest image" magni- 
fication bias is used here. No cut-off due to the finite angular resolution has 
been imposed in this or any of the other figures. 

5.1 Effect of magnification bias 

We also plot the image multiplicities using the total magnification 
bias (which we will abbreviate as b T . In fig. [2] we compare to the 
"second-brightest image" magnification bias (b 2 ). We see that the 
image multiplicities are relatively unchanged, although the lens- 
ing probability for each image configuration is larger by on overall 
factor of 3.5 compared to the b 2 bias (note, however, that this fac- 
tor will depend on the slope of the luminosity function, which we 
have taken as r\ = 2.1). Surprisingly, the ratio of cusps/doubles and 
cusps/quads is slightly larger using the b 2 bias, which is clearly 
shown by the fact that the turning points, when cusps overtake dou- 
bles and quads, are shifted to the left. This is because cusps sys- 
tems nearly always have at least two (and often three) bright im- 
ages, whereas doubles and quads often have one image (or two for 
quads) that is significantly brighter than the others. So using the 
total magnification bias actually boosts the number of quads and 
doubles relative to cusps. 

Any attempt to solve the "quad problem" of radio lenses 
(described in JTJ can then be rephrased as follows: given a cut- 
off at small image separations, does the "turnover point" for 
quads/doubles occur at a sufficiently small image separation such 
that the total number of quads is roughly equal to that of doubles? 
Given the cutoff separation for CLASS, from the graph it is obvious 
that this condition is not satisfied; the majority of lenses still occur 
before 2", where doubles are dominant. It seems reasonable, how- 
ever, that the inclusion of spirals, substructure and external shear 
from lens galaxy environments may push the turnover point to suf- 
ficiently small image separations to solve the quad problem. 

5.2 Comparison to NFW profile and Singular Isothermal 
Ellipsoid 

As a consistency check, it is instructive to compare this distribution 
to that of lensing by galaxies represented by an SIE profile, and 
dark matter halos represented by an NFW profile. These cases are 
plotted in fig. [3] Interestingly, quads and doubles are significantly 
enhanced compared to NFW even at separations as large as 20" . 
It is somewhat surprising that our model agrees so well with 




6 (arcsec) 



Figure 2. Comparison of the "total image" magnification bias (labeled 
bj ) and the "2nd-brightest image " bias (labeled b 2 ). (See j|4.2l for details.) 
We assume a power-law luminosity function, 0(L) ~ IT , The shape of 
all galaxies are taken to be aligned with their host dark matter halo and we 
have fixed galaxy axis ratio to qu = 0.67. 

SIE at low separations. The agreement can be partly understood 
by comparing the elliptici ties of the dark m atter and galaxies. For 
the SIE model we used the lHao et al. ellipticity distribution, 

whereas the galaxy ellipticities in our AC model were directly cor- 
related with their dark matter halos according to q sa i = q + 0.1. 
The q -distribution of the dark matter halos is from Ijing & Sutol 
(2002), and this distribution depends on the mass of the halo. The 
image separation distribution is dominated at low separations by 
10" M G - 10 12 M o systems. We verified that the distribut ion in this 
range roughly coincides with the observed distribution in lHao et"al] 
(2006), though the halo ellipticity distribution has a significantly 
larger tail at low ^-values. Indeed, since naked cusps appear only at 
q < 0.4 for an SIE lens, cusps have negligible lensing probability 
in the SIE model (dP/d® < 10~ 9 ). From the ^-distribution alone, 
therefore, one might expect the AC model to yield slightly more 
quads (and less doubles) than the SIE model. 

However, ellipticity is not the only factor determining the 
quad-to-double ratio. In fact, when compared at identical elliptici- 
ties, our model results in a smaller quad-to-double ratio compared 
to SIE. The reason is that, while the profile of the adiabatically con- 
tracted halo is very close to isothermal, the addition of a Hernquist 
profile actually steepens the slope of the density profile in the vicin- 
ity of the Einstein radius. This enlarges the radial caustic and results 
in more doubles compared to quads and cusps. The net effect is that 
the slightly higher ellipticities in the AC model is nearly balanced 
out by the steeper density profile compared to SIE. In view of these 
considerations, we consider the agreement of our model with SIE 
(at least for doubles and quads) to be somewhat fortuitous. 

5.3 Alignment of projected major axes 

How much does the assumption of alignment vs. misalignment af- 
fect the quad/double and cusp/double ratio? In fig. [4] we plot these 
two cases. (Here, for simplicity we set the galaxy axis ratio equal to 
qi, = 0.67, which is the mean axis ratio in the "Coma model".) The 
quad/double ratio at = 0.3" changes from 0.37 in the aligned 
case, to 0.08 in the unaligned case. The cusp/double ratio changes 
from 6.5 x 10~ 3 in the aligned case to 4.4 x 10~ 4 , over a factor of 
10 difference! We must stress, however, that if galaxies and halos 
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Figure 3. Comparison of the full model to distribution of image multiplicities for an NFW density profile and Singular Isothermal Ellipsoid ( SIE) density 
profile. Going from left to right we have the comparisons for the the total lensing probability, doubles, quads, and cusps. The lensing probability for cusps in 
the SIE model is smaller than 1CT 8 and not shown here. 



are randomly aligned, the image multiplicities will fall in between 
the two extreme cases we have plotted here. The total probabil- 
ity is very similar in either case, changing from 4.8 x 1CT 4 in the 
aligned case to 5.2 x 10~ 4 in the unaligned case (again at = 0.3"). 
The similarity is to be expected since only the overall shape of the 
lensing mass is being changed. The reason for the slight enhance- 
ment in total probability for the unaligned case is due to the mag- 
nification bias; this effect will be further explained later when we 
compare to the lensing probability with spherical mass distributions 

en 

5.4 Effect of adiabatic contraction 

We now investigate what happens if adiabatic contraction is 
"switched off", i.e. if the dark matter halo has an uncontracted bi- 
axial NFW density profile. The results are plotted in fig. [5] Here 
the galaxy and dark matter halo are aligned, and we have fixed the 
galaxy axis ratio at q H = 0.67. The inner slope of the density profile 
is shallower compared to the AC case, resulting in more quads and 
cusps relative to doubles. The point where quads and doubles be- 
come equal has shifted down to = 0.8", so the predicted number 
of quads increases significantly compared to the adiabatically con- 
tracted model. On the other hand, cusps become dominant at 2.5", 
which is inconsistent with the data since at most one naked cusp 
has been observed to date. For comparison, we plot the AC/no-AC 
cases together in fig. [6] Not surprisingly, the lensing probabilities 
in the middle region are markedly higher in the AC case since the 
dark matter halos are more centrally concentrated than their origi- 
nal NFW form. 

We also consider what happens if the axis ratios a/c, b/c are 
increased by a factor of 0.1 following adiabatic contraction. The 
result is qualitatively the same as before, so we will not plot it here; 
however the turnover points are shifted toward higher image sepa- 
rations. Quads dominate over doubles at 3", so their overall fraction 
is predictably reduced. Cusps have shifted from 6" to 10", which 
may be more consistent with the data since, out of the few observed 
lenses with separations greater than 6", none of them have three 
observed images. The quad fraction is obviously less consistent 
with the data in this case. We expect however that numerous fac- 
tors which increase the quad/cusp fraction, such as spiral galaxies 
and galaxy environments, will have a greater effect at galaxy- and 
group-size separations than at cluster-size scales where cusps be- 
come dominant. Thus, models where the overall axis ratios are in- 
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Figure 5. Distribution of image separations when no effect from baryon 
cooling is included. See also Fig.[6]for a comparison of the cases with and 
without adiabatic contraction. The shapes of the projected luminous and 
dark matter densities are aligned and the axis ratio of the luminous compo- 
nent is fixed at qn = 0.67. 

creased by adiabatic contraction may ultimately be consistent with 
the lens data. 

5.5 Comparing to spherical mass distributions 

The difference between the spherical and elliptical case is striking 
as may be ascertained from fig. [7] The dramatic enhancement at 
large separations is a consequence of the triaxiality of the dark mat- 
ter halos; elongation of halos along t he line-of- sigh t create s con- 
siderably larger deflection angles dOguri & KeetorJ J2004h ). The 
enhancement gets larger with increasing separation because larger 
mass halos tend to be more triaxial. The small-separation side does 
not show this effect because we did not start with a triaxial model 
for the central galaxy-rather, the galaxy properties were taken from 
observational data which naturally sees the galaxies in projection. 

Curiously, the lens probability in the spherical case at small 
separations is actually higher than in the elliptical case. This is a 
consequence of the magnification bias. A source placed at the cen- 
ter of a lens with circular symmetry forms an Einstein ring; this 
degeneracy means that sources placed near the center will have 
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Figure 4. Comparison of the distribution of image separations in the cases where the shapes of the luminous and dark matter projected densities are aligned 
and anti-aligned. From left to right we show the distributions for doubles, quads and cusps. The last panel shows all three multiplicities in one plot. We have 
fixed the galaxy (luminous) axis ratio to be qu = 0.67, and the "second brightest image" magnification bias is used here. 
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Figure 6. Comparison of the distribution of image separations when the effect from baryon cooling on dark matter halo (adiabatic contraction) is included 
(dotted) and not included (solid). The shapes of the projected luminous and dark matter densities are aligned and the axis ratio of the luminous component is 
fixed at qu = 0.67. Going from left to right we have the comparisons for the the total lensing probability, doubles, quads, and cusps. 



both of their images highly magnified. A similar lens with ellip- 
ticity, however, does not have this degeneracy. Although sources 
placed inside and near the astroid caustic will be strongly magni- 
fied, a source slightly outside the astroid caustic does not produce 
such high mag nifications since it maps to transition loci rather than 
critical curves iFinchetal.|j2002h . This effect implies that if a lens 
with circular symmetry is deformed into an elliptical lens such that 
the average image separation stays fixed, the magnification bias for 
doubles is drastically reduced while the high magnification bias for 
quads may or may not compensate for this reduction. Whether the 
total bias cross section is increased or reduced by ellipticity de- 
pends on the type of bias used, and on the density profile of the 
lens. For our lens model, if the "second brightest image" bias is 
used (as in fig. [7]l the total cross section will be lower than in the 
spherical case, whereas if the total magnification bias is used it will 
be roughly the same. Hence, depending on which type of bias is 
appropriate for a set of observed lens systems, the assumption of 
spherical galaxies may result in a slight overestimate of the total 
lensing probability. 
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Figure 7. Comparison of the distribution of image separations in models 
with spherical and elliptical projected mass profiles. For the spherical case, 
we separately also include scatter in the mass-luminosity and concentration- 
mass relations. We plot the total lensing probability in the triaxial case 
(same distribution as fig. [TJ for comparison - see curve labeled "ellipti- 
cal". The "second-brightest image" magnification bias is used for all the 
calculations in this plot. 
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5.6 Scatter in mass-luminosity and concentration-mass 
relation 

In calculating all the results presented thus far, we ignored the ef- 
fect of scatter in the mass-luminosity relation and in halo concen- 
trations. The question naturally arises as to how scatter effects the 
lensing probabilities. To investigate this we calculate the total lens- 
ing probability assuming spherical galaxies and halos. The result is 
fig. [7] with the total lensing probability in the elliptical case (from 
fig- IB being plotted for co mparison. For the sp herical concentra- 
tions of halos we used the iBulIock et al.l ( 120011) model, assigning 
log-normal scatter around the median with standard deviation cr = 
0.3. We also assumed log-normal scatter i n mass-luminosity with a 
standa rd deviation S = 0.25, as derived in lCoorav & Milosavlievid 
d2005h . Scatter in mass-luminosity has the effect of increasing the 
total lens probability by «11%. Scatter in concentration increases 
the contribution of the dark matter to the lensing probability, and 
has a noticeable effect at separations as small as 2" . While includ- 
ing scatter in these relations will not have a profound effect on im- 
age multiplicities, the concentration scatter will probably increase 
the quad fraction slightly in the middle region (as 2" -10") since the 
dark matter halos produce higher quad and cusp fractions than the 
galaxies. 



6 DISCUSSION 

We have shown that the distribution of strong lensing image mul- 
tiplicities offers new possibilities for constraining correlations be- 
tween galaxies and their host dark matter halos — in particular, adi- 
abatic contraction and the degree of axial alignment. Since varying 
the amount of alignment makes a considerable difference in the 
quad/double ratio, and given that the observed ratio is so high, one 
might expect galaxies and h alos to be qu i te closely aligned. Simu- 
lations seem to support this iB ailin et al.l j2005l) .iKazantzidi s et al.l 
and indeed there is some observational evidence from lens- 
ing: iKochanekl (2002) used mass profile modeling of 20 lenses 
from the CfA-Arizona Space Telescope Lens Survey (CASTLES) 
to show that the mass and light distributions were aligned within 
(Ad 2 ) 112 < 10° where 6 is the projected alignment angle. 

The distribution of image multiplicities provides a simple 
method, far less time-consuming than individual lens modeling, 
that can be applied to a large statistical sample of lenses. This 
method also has the advantange that the data for each individual 
lens does not have to be very high quality. The disadvantage is that 
we have to include all correlations and enviromental effects that are 
either constrained by other data or directly obtained from theory. 

One important caveat to keep in mind is that alignment of 
the major/minor axes of galaxy and halo does not necessarily im- 
ply that the projected axes will be aligned. The projected axes can 
also be misaligned if the galaxy and halo have differing triaxialities 
( lKeetonetai] ( ll997l) ). analogous to the phenomenon of "isophote 
twist". Our method constrains the alignment of the projected axes. 
If the projected axes are shown to be closely aligned, then we can 
draw the strong conclusion that the major/minor axes of galaxy and 
halo are aligned, and that the triaxial galaxy and halo shapes are 
similar on average. 

There are several ways in which our model can be improved. 
First, a fully triaxial model for adiabatic contraction from sim- 
ulations is required . We used analytic fitting functions from the 
Gned in et al ] 120041) model of adiabatic contraction that was de- 
rived from halos with spherically averaged density profiles, extend- 
ing it to triaxial halos by making the replacement r — » R where R 



is the triaxial radius. But this is clearly only an approximation that 
may break down at high ellipticity. The predicted lensing probabil- 
ity for cusps may differ significantly when a more accurate model 
of adiabatic contraction is used. 

In addition, we have neglected the contribution from spiral 
galaxies. This is partly justified by the fact that at least 75% of 
galaxy lensing is thought to be due to ellipticals ( Moll er et al.l 
(2007)). However, since disk galaxies viewed edge-on can have ex- 
tremely high projected ellipticities, their effect on image multiplic- 
ities may be substantial. We have also not included scatter in the 
mass-luminosity relation and in halo concentration, the effects of 
which were discussed in the previous section. Although the image 
multiplicities are not drastically effected by the scatter, it is nev- 
ertheless an important effect and has an appreciable impact on the 
total lensing probability. 

We have also ignored the effect of external shear from lens 
galaxy environments, which can incre ase lensing probabilities at 
intermediate image separations, > 3" dOguril 420051) ). Another ef- 
fect is the misalignment of the center-of-mass of galaxies with 
that of their host halos, which will s lightly enhance the cross sec- 
tion for quads ( Ouadri et al. I J2003k Perhaps most importantly, 
we have not included su bstructure in our calculations. Work by 
ICohn & Kochanekl te004h suggest that satellites can roughly dou- 
ble the expected quads-to-doubles ratio, thus dramatically altering 
the predicted ratio from our model. They also give rise to cross sec- 
tions for more exotic image configurations, with five or more visi- 
ble images, whose lens probabilities may be significant. Including 
all these effects will be essential for deducing rigorous constraints 
on the profile, shape and correlations of lensing masses from future 
data. 
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